Unified dark energy and dark matter from a scalar field different from quintessence 
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We explore unification of dark matter and dark energy in a theory containing a scalar field of non- 
Lagrangian type, obtained by direct insertion of a kinetic term into the energy-momentum tensor. 
This scalar is different from quintessence, having an equation of state between —I and and a zero 
sound speed in its rest frame. We solve the equations of motion for an exponential potential via a 
rewriting as an autonomous system, and demonstrate the observational viability of the scenario, for 
sufficiently small exponential potential parameter A, by comparison to a compilation of kinematical 
cosmological data. 
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I. INTRODUCTION 

Dark energy remains a fundamental mystery, both in 
terms of its unexpectedly low but non-zero value and 
because of the apparent coincidence of its present den- 
sity being approximately that of other components. At- 
tempts to address the coincidence problem have been of 
two types. One is to invoke the anthropic principle, per- 
haps at its most persuasive when coupled with the con- 
cept of the string landscape [l[ . The second is to permit 
the dark energy to be a dynamical entity, and hope to 
exploit solutions of scaling or tracking type to remove de- 
pendence on initial conditions (for a review see Ref. (H). 
It is however fair to say that no compelling scenario of the 
second type has been found that is compatible with the 
tight present observational constraints on the equation of 
state parameter w 

In this paper we study the consequences of model- 
ing the dark energy using a scalar field that is of non- 
Lagrangian type. The principle that fundamental physics 
should derive from a Lagrangian description is a deep- 
seated part of modern physics, as powerfully argued for 
instance by Durrer and Maartens [4J . It is a measure of 
the difficulty of the dark energy problem that there have 
been several papers that have abandoned this principle, 
for instance modeling the dark energy as a phenomeno- 
logical fluid which exhibits a particular scaling with the 
scale factor [f| or Hubble parameter [f| , or even allowing 
a cosmological constant with an explicit dependence on 
time Q- Our proposal too is of this general type; we are 
closer to traditional quintessence modeling in adopting a 
scalar field description, but consider a scalar field that 
does not emerge from a Lagrangian. 

Dropping the Lagrangian assumption is a major step, 
and in taking such a step one wishes to be sure that 
there is significant payback. Our model offers one such 
reward — it permits a unified description of dark energy 
and dark matter as due to the single field we consider. 
While our proposal is a speculative one, this opportunity 
is significant enough to merit study. 



Our proposal is not of course the first to seek to unify 
dark energy and dark matter into a single material. Dis- 
counting those where the dark energy arises from a con- 
stant term in the action, some examples are as follows. 
Ref. Q proposed a tachyon-type scalar-field Lagrangian, 
in which the scalar fluid can be broken up into dark mat- 
ter and dark energy components. K-essence unification of 
dark matter and dark energy has been studied in Ref. M. 
Staying instead with the canonical Lagrangian, Ref. [10( 
introduced a complex scalar field with a mixed potential 
made of quadratic and exponential terms, which then 
mimic dark matter and dark energy respectively. Alter- 
native strands with similar goals are study of the gener- 
alized Chaplygin gas [ll| and of barotropic fluid models 

m. 



II. MODIFYING THE EINSTEIN EQUATIONS 



The usual Einstein equations are given by 



(1) 



where G^, T^ v and Ag^ u are the Einstein tensor, 
energy-momentum tensor and the cosmological constant 
term, respectively, is the metric tensor, A is the cos- 
mological constant and k 2 = 8ttG. A is the simplest ver- 
sion of dark energy, being time independent and isotrop- 
ically and homogeneously distributed in space. It suffers 
from the coincidence problem, and to address this we 
wish to allow the cosmological constant to evolve. 

A simple idea, as adopted in quintessence models, is 
to allow A to be a function of some scalar field <j>. One 
cannot however allow this dependence to be on <f> alone; 
the Bianchi identity 

VG MV = , (2) 
and the law of energy-momentum conservation, 

V%, - , (3) 
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would force 



III. DYNAMICAL EQUATIONS 



V^ 1 [A(«^]=0 



(4) 



requiring A(</>) to still be constant. 

To remedy this it is necessary to incorporate a dynam- 
ical term, depending on V M </>, into the equations. For 
quintessence this is done by including a canonical kinetic 
term in the Lagrangian; A((f>) then becomes the scalar 
field potential and the total dark energy density includes 
both potential and kinetic terms. Here we propose the 
simplest possible alternative, which is the direct insertion 
of a kinetic term into the energy-momentum tensor: 



Gnu — K 



+ A {4>) g^ - iv^V, 



(5) 



Now A(<j>) is not necessarily a constant. The equation of 
motion for the scalar field is then given by 



dA V^V"( 



(V M V t 



= o, 



(0) 



This equation follows directly from the Einstein equa- 
tions (plus the assumption that for the other components 
T^ v remains separately conserved) , but this form is more 
convenient. 

At first glance this scalar field looks very much like 
a quintessence field. But in fact it is very different 
from quintessence, and even different from K-essence (l3j 
where the Lagrangian is written as a general function of <p 
and V^4>V ^4>. Indeed, it has no Lagrangian formulation 
in the framework of K-essence theory. 

We can present the proof as follows. In general, the 
Lagrangian of K-essence is given by an arbitrary function 



where 6 is a scalar field and 



X 



1 



(7) 



(8) 



Since we choose the signature of (— 1,+1,+1,+1), we 
always have X > 0. Varying this Lagrangian with respect 
to the metric we obtain the energy-momentum tensor in 
the form 



(9) 



where _£? x denotes partial derivative of the Lagrangian 
with respect to X. By identifying it with the energy- 
momentum tensor of the scalar 

T^u = --V^V^ + A (0) , (10) 

we find the corresponding Lagrangian does not ex- 
ist. This demonstration is limited to Lagrangians 
of K-essence form, but there is no reason to 
think that a more general Lagrangian, such as 
I£ (X, (f>, R^V^fiVcj), ■ ■ ■ ) could lead to our equations 
while retaining the Einstein form of gravity. 



The equations that govern the evolution of the 
spatially-flat Friedmann-Robertson- Walker Universe are 



3H 2 = n 2 



2H + 3H 2 



1 ;< 



-A(0) + ^J , 



^ + 3tf (jH+Pi)=0, 



(11) 



(12) 



where H = a/a is the Hubble parameter, pi and pi are 
the energy density and pressure of i-th matter compo- 
nent (namely, relativistic matter, baryonic matter, and 
so on), a is the scale factor and dot denotes derivative 
with respect to the physical time t. We have set a^ = 1 
for the present universe. 

From the Einstein equations above we obtain the den- 
sity and pressure of our scalar 



Psca = ~<j> 2 + A (0) , 
p sca = -A (<p) . 



(13) 



These can be contrasted with the equivalents for 
quintessence with the same potential 



Pqui = -<t> 2 + A ((/)) , 
Pqui = - A 0) • 



(14) 



From the expressions of density and pressure, we know 
quintessence has the equation of state — 1 < w qu i < 1 
for A > 0, while the scalar has — 1 < w sca < 0. From 
the conservation equation (|12p we then know that the 
density of quintessence scales in the range a~ 6 to a , 
while for the scalar the range is restricted to a~ 3 and 
a . 1 This property suggests that the scalar may play the 
role of both dark matter (scaling approximately as a -3 ) 
and dark energy (scaling approximately as a ). 

From the expressions for the density and pressure we 
can further derive the sound speed in the rest-frame of 
the scalar fields, and find that they are different as well: 

_ dp/dX 



d P /dX 



. 2 dp/dX _ 

Cq U ; „ , 1 



dp/dX 



(15) 



where X = </) 2 /2. This is the case for any potential A(^>). 
The vanishing of the sound speed allows our scalar field 
to cluster gravitationally more easily than quintessence. 



1 This is the same range accessible to the simplest DBI tachyon 
model but our equation of motion differs from that case. 
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Since this is crucial in order to match the cosmic mi- 
crowave background (CMB) data, we derive the behavior 
of the perturbations explicitly later on, confirming this 
result. 

The equation of motion for <f> can be derived from 
Eqs. (HI]) or Eq. © as 



(16) 



Compared with the equivalent equation for quintessence, 
with the same potential, there is a significant difference: 
the friction term in the equation of motion for the scalar 
is only half that of the quintessence. So with increasing 
redshift, the densities of the scalar field will increase more 
slowly than quintessence. In fact, if the potentials are 
constant or sufficiently flat such that dA/d(f> ~ , in a 
kinetic-dominated regime we will have 



Pqui cx a , 



for quintessence and 



Psca oc a , 



(17) 



(18) 



for the scalar. The latter is exactly that of cold dark 
matter. 



IV. THE EVOLUTION OF THE 
HOMOGENEOUS SCALAR FIELD 

In order to compute the evolution of the scalar field and 
to check whether it is compatible with current data sets, 
we need to specify a form for the potential A(4>). The 
simplest choice is just a constant potential, A = const in 
Eqs. (fTTj) . Then from Eq. (|16[) we have the density 



A 



PdO 
„3 



(19) 



where pdo is a constant which can be interpretted as 
the present dark matter density. This is exactly the 
ACDM model, with A playing the role of dark energy 
and V^V u cf>/ 2 the role of dark matter. 

However, taking the potential to be constant is effec- 
tively reintroducing a pure cosmological constant (c.f. 
the kinetic K-essence model of Scherrer [9j), and hence 
does not represent a significant step forward in under- 
standing the nature of dark energy, although our proposal 
has a novel nature for the dark matter. We therefore 
choose a more general form of the potential, 



M4) = V t 



-K\(p 



(20) 



which we will use throughout the remainder of the paper. 
Here Vq and A are two constants. Without loss of gen- 
erality, we assume A > 0. In the limit A — > we recover 
the constant potential case and therefore the model is 
continuously connected with ACDM, at least where the 
background evolution is concerned. 



A. Autonomous system of equations 

To study the evolution of the field, we set up an au- 
tonomous system. The main equations are given by 



3H 2 = n 2 



\^ 2 + A (</>) + p t + Ph 



2H + 3H 2 = -k 2 



-A {<t>) + ~ Pl 



-H(f> 



dA 



= 



(21) 
(22) 
(23) 



where p r and pb are the density of radiation and baryonic 
matter, respectively. They have a constant equation of 
state equal to 1/3 and 0, respectively. Following Ref. [ll|, 
we introduce the following dimensionless quantities 



K(p 



V6H ' 



y ; 



ka/A 
V3H ' 



(24) 



Here x 2 and y 2 represent the density parameters of the 
kinetic and potential terms respectively. We expect in- 
teresting cases to have the scalar field rolling down the 
slope of the potential, so since we have assumed A > 0, 
we should have x > 0. Then the above equations can be 
written in the following autonomous form 



dx 
dN 



dy 
dN 



dilh 
~dN 



3 

+ 2 X 



V6 



1 



5<' 



H — V 
-30 b 



i-y 2 



x 2 -y 2 



-(i-x 2 - y 2 - n h ) 



together with a constraint equation 



y 



tt h + f2 r = 1 



(25) 

,(26) 
(27) 

(28) 



Here N = In a. The equation of state and the fraction 
of the energy density f2^, for the scalar field are 



y 



ZH 2 



1 i 2 

x +y 



(29) 
(30) 



B. Observational requirements 

What constraints does a model seeking to explain both 
dark matter and dark energy have to satisfy? Normally, 
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Name 


X 


y 




Existence 


Stability 






(a) 











All A 


Unstable node 







(b) 


X 





Vl-x 2 


All A 


Saddle line segment for < x < 1 


x 2 





(c) 


^A 







A 2 <| 


Stable node for A < ^ 


1 


-1 + §A 2 



TABLE I: The properties of the critical points for the exponential potential given by A = Vbe 



observational constraints are imposed under the assump- 
tion of separate dark matter and dark energy compo- 
nents, but due to the dark degeneracy, first described 
by Hu and Eisenstein [l6[ and then further explored in 
Refs. [13, EM Gil, gravitational probes alone are unable 
to give a unique decomposition and can only impose con- 
straints on the total dark sector. In Ref. [2fj] we re- 
cently derived the constraints on a combined dark sector 
fluid from current kinematical observations in a model- 
independent way. These showed that the total dark sec- 
tor equation of state must start at or near the cold dark 
matter value w = 0, and then evolve to become nega- 
tive by the present following a particular profile. The 
standard cosmological model, e.g. as in Komatsu et al. 
Q , predicts a present total dark sector equation of state 
of about —0.78 (the weighted mean of the dark matter 
and cosmological constant contributions), but in fact this 
value is only weakly constrained [13, HH . The behavior 
is more tightly constrained at higher rcdshifts, where the 
actual observational data lies. In addition, a success- 
ful unified dark sector model must reproduce the present 
dark sector density O^ark = 0.96. Our aim will be to test 
whether our model can achieve this. 



C. Fixed points and phase portraits 

In Table I, we present the properties of the three fixed 
points for the exponential potential. The point (a) corre- 
sponds to the radiation-dominated epoch and this point 
is unstable. The line segment (b) corresponds to a scalar 
plus baryon-dominated epoch and it is a saddle line seg- 
ment. In this epoch, the scalar field behaves as dust 
matter which has the equation of state w — 0. The point 
(c) corresponds to a scalar-dominated epoch. Point (c) 
is stable and thus an attractor. In this epoch, the scalar 
has an equation of state w < — 1/3 if A < 1, and so the 
Universe accelerates in this epoch. 

Viable scenarios start at high redshift near the unsta- 
ble radiation fixed point (a). This is necessary since the 
distance from the origin corresponds to the relative en- 
ergy density in the scalar field. As that energy density, 
like the one in matter, decreases slower than the radia- 
tion energy density, we need to start close to x — y = 0, 
analogous to the 'thawing' regime of quintessence. In 
order to follow the usual evolution of the Universe, the 
field should then move to the effective matter-dominated 
saddle line. This happens if we start with y <C 1 as 
dy/dN oc y w in this limit. The trajectories then move 




FIG. 1: The phase plane for the best-fit A found below, 
approximately 0.1. The point (0, 0) corresponds to the 
radiation-dominated epoch. The point (0, 0) is unstable and 
the point (0.0082, 0.9999) is stable and thus an attractor. 
The line segment(a;, 0) is a saddle line. The initial con- 
ditions best matching observations lead to the red/dotted 
trajectory, while the other trajectories have different ini- 
tial conditions. The outer thin solid line corresponds to 
x 2 + y 2 = 1 — flbo — firo, giving the corrrect present dark 
sector energy density, while the angle B gives the required 
equation of state w$ = — sin 2 B ~ —0.78. So the present-day 
Universe must lie in the vicinity of the point A which is the 
intersection between the two thin solid lines, which the dotted 
line indeed passes through. 



across, staying close to the y = line, before turning up. 
Figure [1] shows phase portraits for A = 0.01 with various 
initial conditions. The trajectories are all confined inside 



the circle given by x 2 + y 2 = 1 due to the constraint 
equation Eq. ((28)) . 

The present epoch can be identified through the re- 
quirement that flfio = x 2 + y 2 ~ 0.96, which corresponds 
to a circle just inside the limiting circle x 2 + y 2 = 1. 
Since all trajectories evolve towards the single attractor 
(c) which lies on x 2 + y 2 = 1, all viable trajectories will 
cross that line eventually. In addition, the total equation 
of state parameter of the scalar should be of the order 
of wo — —0.78 today. That condition can be graphi- 
cally represented by a straight radial line with an an- 
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FIG. 2: The evolution of density fractions for radiation 
(magenta/dot-dashed), baryons (blue/dashed), kinetic term 
x 2 (black/solid), and potential term y 2 (red/dotted), for the 
best-fit model found in subsection II V Dl 

gle B ~ arcsin(V0.78) with respect to the a;-axis since 
W0 = — sin 2 (B). The good models then need to cross the 
circle of today's at the intersection with this line. We 
will examine the observational constraints in more detail 
in the following subsection. 

In Fig. [2j we plot the evolution of density fractions 
for radiation, baryon matter, kinetic term and potential 
term, for best-fit model parameters we determine below. 
This shows that the scalar field can mimic the cold dark 
matter and dark energy very well. 



D. Constraints from current data 

We now impose detailed observational constraints on 
our model to establish its viability. We follow a method 
very similar to that outlined in Ref. [2(| to constrain the 
model, applying a Markov Chain Monte Carlo (MCMC) 
approach to compute the parameter posterior probabil- 
ities. We assume a flat universe, and fix the radiation 
density today from the CMB temperature. The evolu- 
tion of the scalar field is defined in terms of the potential 
parameter A and the value of the equation of state today 
wo , and then integrated backwards to finds its evolution 
at earlier time. We assume a uniform prior on the equa- 
tion of state parameter of — 1 < wq < and a log prior 
on A such that —4 < log 10 (A) < 1. We also include the 
baryon density flbh 2 and the Hubble parameter today 
Hq as free parameters, and marginalize over them. 

We use a fairly typical compilation of kinematical data. 
Standard candle data comes from supernova type la lu- 
minosity distances, for which we use the cut Union su- 
pernova sample [22j (with systematic errors included), 
and standard ruler data comes from the angular posi- 
tions of the CMB [23| an d Baryon Acoustic Oscillation 
peaks 24]. Note that Ref. [23] gives constraints on the 
scaled distance to recombination R and the angular scale 




-0.8 -0.75 -0.7 -0.65 -0.6 -4 -3 -2 -1 1 

w„ log 10 « 



FIG. 3: The 1-d and 2-d probability distribution for the equa- 
tion of state today, wo, and the potential parameter, log 10 (A), 

of the sound horizon Z a . These are defined to be 

R^J^t(zcmb), l^^p^l- (31) 

Since R is scaled by the physical matter density, and 
so makes assumptions about the separability of the dark 
matter and dark energy, we ignore it in this work. We 
use only the constraints on Z a , as well as those on Vl^ti 2 
and the correlations between the two. We also include 
the SHOES [25[ measurement of the Hubble parameter 
today, H = 74.2 ± 3.6 kms^Mpc -1 . 

We find that the model is a good fit to the data. The 
best-fit parameters have a x 2 = 312.1, which is almost 
equivalent to the best fit of the LCDM model, x 2 = 311.9 
(though the scalar field model has one extra parameter) . 
The equation of state today lies in the range —0.82 < 
wo < —0.57 at 95% confidence. The 95% upper limit 
on the potential parameter is A < 0.20. The probability 
distributions for these two parameters are shown in Fig. [3] 
and some sample w(a) curves in Fig. 

V. STRUCTURE FORMATION 

From the analysis in Ref. [lj| we know that the model 
will fit the CMB data if the rest-frame sound speed is 
indeed zero. Because this is such an important condition 
on the model, we here derive the sound speed directly 
from the perturbation equations. We will discuss the be- 
havior of the perturbations further in the Appendix. For 
this purpose, we work in Newtonian gauge. In the ab- 
sence of anisotropic stress, and for scalar perturbations, 
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we get the perturbed counterpart 



• • .2 ciA 

= O £0 - $00 + -j-S(j) , 



dA 



2a 



00(50 = p^o^ 



02? 



- (i, x) 8\ 



(38) 



FIG. 4: The evolution of the equation of state for the scalar, 
shown for 8 models drawn from the Markov chain. It behaves 
as the cold dark matter at higher redshifts and dark energy 
for the lower redshifts. For comparison the evolution of the 
total w in the ACDM model is shown as the thicker black line, 
which runs more or less centrally through the set. 



the perturbed Friedmann-Robertson- Walker metric can 
be written in the form 



Since we are working in linear perturbation theory, it is 
convenient to transform the equations from real space to 
Fourier space since each Fourier mode evolves indepen- 
dently. We will also suppress the O-subscripts for the 
homogeneous and isotropic quantities from now on. 

It is well known that both the adiabatic sound speed 
and the rest frame sound speed (the sound speed for the 
fluid in its rest frame) play a very important role in the 
discussion of structure formation theory. Here we work 
out the two quantities explicitly. The adiabatic sound 
speed squared is defined through [26] 



2 dA 



p<j> 3iJ0 



(39) 



ds 2 



(1 + 2$) dt 2 + a (t) 2 (1 - 2$) dx i dx i , (32) 



where $ is the gauge-invariant Newtonian potential. The 
potential characterizes the metric perturbations. 

For the dark matter and dark energy dominated Uni- 
verse, we can safely neglect the effect of radiation. We 
assume baryonic matter as a perfect fluid which has the 
energy-momentum tensor 



The rest frame sound speed squared c 2 of the scalar is 
related to the pressure perturbation in the Newtonian 
gauge through 



S Pcj> = clSptj, + 



3aH 



(40) 



Expressing this equation with perturbation quantities, 
following the procedure in Ref. [26|, we find 



Tftu = {fib + Pb) U^U V + Phg^y 



(33) 



• , dA ■ \ dA 

2 $ + —50 + 3^0,50 - — £0 



where is the four-velocity of the fluid. Perturbations 

in the energy density p h , pressure p h and four- velocity Therefore, we can conclude immediately that 



can be written as 



Pb (t, x) 
Pb (t, x) 
tV (t, x) 



Pob + 5pb (*, x) , 
Pob + 5p h (t, x) , 
(°)u M + <5u M (t,.x) , 



cf =0. 



(41) 



(42) 



(34) 



where ^\ 



(—1, 0, 0, 0) and pob{t), Pob{t) are the 



This is consistent with the discussion in Eq. (15). The 
vanishing sound speed in the rest frame of the scalar al- 
lows it to play the role of cold dark matter. 

We have the perturbation for the pressure as follows 



homogeneous and isotropic energy density and pressure. 
So we obtain 



,3aH 



k 2 



-pV. 



(43) 



6T ° = Sp h (t, x) , 

STq = (pob + Pob) Su l (t,x) , 

STl = -6 Pb (i, x) 51 . 

For the scalar field, we define the perturbation as 

(t, x) = 0o (t) + <50 (t, x) . 
From the energy-momentum tensor 
1. 



rV M 0V 1/ + A(0) ff(tl „ , 



(35) 



(36) 



(37) 



Both CDM and a cosmological constant have 5p = 0. 
How large is the contribution to dp which arises from the 
gauge transformation to Newtonian gauge due to c 2 ^ 0? 
The adiabatic sound speed squared can be written as 



W<f, 



3H (1 + w ) 



(44) 



We note that the adiabatic sound speed is determined by 
the homogeneous quantities. To have a picture of the adi- 
abatic sound speed, we should resort to the background 
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FIG. 5: The adiabatic sound speed squared (solid line) 
and the equation of state (circled line) for the scalar, showing 
c:; ~ Wj,. At redshifts greater than 6, they are both vanishing. 

equations, given by 

H = ~ (Pb + ^ 2 ) , (45) 

+ |^ + |=Q. (46) 

Using our best-fit parameters, in Fig. [SJ we plot the equa- 
tion of state w<j> and the adiabatic sound speed squared 
c 2 for the scalar. It shows that ~ w$. At the redshifts 
greater than 6, both the equation of state w$ and the 
adiabatic sound speed c a are nearly zero and unimpor- 
tant for structure formation. This point is essential for 
the scalar to play the role of cold dark matter. In conclu- 
sion, we have c a ~ 0, ~ at redshifts greater than 
about 6. Thus we also have Sp ~ at redshifts greater 
than 6, from Eqs. (gSJ and ijlijl. 

The Appendix further explores the properties of the 
structure formation equations. 

VI. CONCLUSION 

We have investigated in a cosmological context the be- 
havior of a scalar field with non-standard kinetic term in 
the Einstein equations. So far we lack a Langrangian de- 
scription for this scalar, and we have shown that it is not 
possible to build one in the framework of K-essence fields. 
We have not been able to exclude that a more general 
Lagrangian, such as J?(X, (j), i? M ^V M 0V ly 0, • • •) could 
mimic our equations while retaining the Einstein form of 
gravity, or that the Lagrangian looked for could exist in 
the framework of Kaluza-Klein theories, but equally we 
have no reason to think it will. 

On the plus side, we find that this scalar has some in- 
teresting properties. In the first place, it has an equation 
of state between w = — 1 and w = 0. This is differ- 
ent from the quintessence field which has the equation 
of state between w — — 1 and w = +1. Hence the scalar 



field can behave as pressureless matter in the matter- 
or radiation-dominated epochs, later evolving to take on 
dark energy properties as well. A degree of fine-tuning 
is needed in the initial conditions in order to ensure that 
the scalar field only dominates in the latter stages of the 
evolution, which is of the same form as that invoked in 
thawing quintessence models. 

Secondly, the rest frame sound speed of the scalar is 
zero. This is different from quintessence for which the 
sound speed is equal to the speed of light. Although a 
DBI scalar field has the same range of equation of state as 
our scalar, its rest frame sound speed is nonvanishing. As 
is known, a vanishing sound speed is sufficient for a scalar 
field to play the role of cold dark matter in the process 
of structure formation, and the sound speed should not 
be too big if gravitational collapse is to match current 
observations pjl]. Thanks to the vanishing sound speed 
of our scalar, we find that it behaves exactly as cold dark 
matter in the process of structure formation. 

To conclude, we have invoked a scenario in which a 
non-Lagrangian scalar field is able to play the roles of 
both dark matter and dark energy. With sufficient tuning 
of initial conditions, we have shown that a satisfactory 
evolution can be arranged, with present data constrain- 
ing an exponential potential to have an exponent of 0.2 
or less, though the model does not significantly improve 
on the fit of the ACDM model. Our model would be 
falsified in the event of direct detection of conventional 
dark matter particles. Whether it can be motivated in 
terms of fundamental theories of physics, perhaps in an 
effective theory formulation, remains to be seen. 
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APPENDIX A: EVOLUTION OF THE 
PERTURBATIONS 

In order to make numerical calculations, we should 
rewrite the equations in the dimensionless form. To this 
end, we use the variable N 

d d d 2 d 2 d 

di =H dN' W = H2 dN 2 ~ + HH 'dN- (A1) 



8 



Here prime denotes the derivative with respect to N. So 
the background equations become 



HH' 



H 2 4> 



Pb 



2 v 



- (3H 2 + 2HH') 4>' 



(A2) 

dA , , . 

— = . (A3) 



Define 



; H o kVo 



_iffw (A ., 



where pbOj ar e the energy density of baryon mat- 
ter and the Hubble constant in the present-day Universe. 
Using these new variables, we can rewrite the main equa- 
tions in the dimensionless form 



(A5) 



ft V + 77 (3/i + 2ft') tf/ - 3-fi A e" KA0 = . (A6) 

2 K 



for baryon matter [27|, and 



Scb = 



= - (1 + W(f> ) ( ^ - :-!<!> 



-3H 



w<t> } 04, , 



= -#(1-3^)0^ 



1 + W0 



,,/fc 2 , kK 

1 + a a 



for the scalar. Here 



$P4> 



2$ 



dA; 



-A 



¥ 2 



A 



(A14) 



(A15) 



(A16) 
(A17) 



The perturbation equation for the scalar is given by 



We can absorb k into </>. Then the main equations are 
simplified to be 



hh' = -\£l hQ e-* N -\h?<t>'\ (A7) 
ftV + \ (3ft + 2ft') 0' - 3Aft A e~ A0 = . (A8) 



■S**+^»£-!**+$W-0.(AU) 



We find it is convenient to consider the following equa- 
tions 



The perturbed Einstein equations are [27 



+ + — * = -y (PbSb + P^), (A9) 

K 2 

k 2 H$ + fc 2 $ = —a [(p h + Pb ) b + fa + P4> ) 64} , (A10) 

2 \ _ k 2 f dA , 



$ + + ( 2i? + 3i? 2 ) $ = y ( ) , (All) 



where <J b ,0 = 3pb,4>/ Pb,<f> is the density contrast for baryon 
matter and scalar, respectively, and 6^,0 = ik ■ Vb,(j> rep- 
resents the divergence of velocity for baryon matter and 
scalar, respectively. We note that 5^, is different from the 
perturbation of the scalar, 8(j>. 

On the other hand, the energy conservation equation 
(which includes the continuity and Euler equations) holds 
for the baryon matter and the scalar field, respectively. 
So we obtain [III 



$ + 4H $ +12H + 3H 2 $ 



3 • k 2 r dA 5 

S <t , = -(l + w <f ,) 

- 3H (I; " w *) S * ' 



It'll, 



1 + 



+ 



k 2 k 2 
1 + w<f, a a 



Using the definition of 



H , (A19) 



(A20) 



(A21) 



(A22) 



(5 b = -—+3$, 
a 

fc 2 $ 

6\j = —HO^ H 



(A12) 
(A13) 



-no -no 



(A23) 



and the variable N, we can rewrite the above equations 
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in the dimensionless form 

,2*" i I AU2 



h 2 <$> + (4ft 2 + hti) $' + (2hti + 3ft 2 ) <$> 



3 



4' + \ {3h 2 + 2hti) 54>' + ^Scj)- 
-eXttxe-^® + 3fl\X 2 e^ x ^S(j) = 



».; = -(i+^)( ^-3*' 



e 



/ = - (1 - 3w0) 6^ - - 
1 + aft aft 



W0 



(A24) 



-ft^$ 



(A25) 



(A26) 



(A27) 



To be consistent with the discussions of the background 
equations, we have rescaled <j) by 4>/k as was done earlier. 
Correspondingly, we have here 



Now we have <E>, S(j>, 5$, and 0^, totalling four pertur- 
bation variables, and four differential equations, namely 
Eqs. (|A24|) - (|A27|I . Thus the system of equations is closed. 



At redshifts greater than 6, we have w^, ~ 0, c 2 ~ 0, 
and Sp^/Spcf, ~ 0. So the perturbation equations simplify 
to 



ftV 



<v = 



(4ft 2 + ftft') $' + (2ftft' + 3ft 2 ) $ = ,(A30) 



ah 



«2 



(A31) 
(A32) 



■fft 2 ^^ 



^ 3(1+^^,) 



ft 2 ^'50' - ftV 2 $ - 3Afi A e- A ^(S0 



. (A28) 
(A29) 



These are none other than the perturbation equations 
for cold dark matter in ACDM model. Therefore, the 
scalar really behaves as cold dark matter in the process 
of cosmic structure formation. 
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